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1. INTRODUCTION 

1.1. The CLEO III Detector 

The CLEO III detector was designed to study 
decays of b and c quarks, t leptons and T mesons 
produced in e + e~ collisions near 10 GeV center- 
of-mass energy. The new detector is an upgraded 
version of CLEO II @. It contains a new four- 
layer silicon strip vertex detector, a new wire 
drift chamber and a particle identification sys- 
tem based on the detection of Cherenkov ring 
images. Information about CLEO III is available 
elsewhere [^,||. 

CLEO II produced many physics results, but 
was hampered by its limited charged-hadron iden- 
tification capabilities. Design choices for par- 
ticle identification were limited by radial space 
and the necessity of minimizing the material in 
front of the Csl crystal calorimeter. The Csl im- 
posed a hard outer radial limit and the desire 
for maintaining excellent charged particle track- 
ing imposed a lower limit, since at high momen- 
tum the error in momentum is proportional to the 
square of the track length. The particle identifi- 
cation system was allocated only 20 cm of radial 
space, and this limited the technology choices. 
We were also allowed a total material thickness 
corresponding to only 12% of a radiation length. 



2. DETECTOR DESCRIPTION 
2.1. Detector Elements 




Figure 1. Outline of CLEO III RICH design. 

The severe radial spatial requirement forces a 
thin, few cm detector for Cherenkov photons and 
a thin radiator. Otherwise the photons have lit- 
tle distance to travel and it becomes very difficult 
to precisely measure the photon angles. In fact, 
the only thin photon-detectors possible in our sit- 
uation were wire chamber based, cither Csl or 
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a mixture of triethylamine (TEA) and methane. 
Use of Csl would have allowed us to use a liquid 
freon radiator with quartz windows in the system 
using the optical wavelength region from about 
160-200 nm. However, at the time of decision, 
the use of Csl was far from proven and, in any 
case, would have imposed severe constraints on 
the construction process which would have been 
both difficult and expensive. Thus we chose TEA 
+ CH4 and used Cherenkov photons between 135- 
165 nm generated in a 1 cm thick LiF crystal and 
used CaF 2 windows on our wire chambers (LiF 
windows were used on 10% of the chambers) Q. 

Details of the design of the CLEO III RICH 
have been discussed before g. Here we briefly 
review the main elements. Cherenkov photons 
are produced in a LiF radiator. The photons 
then enter a free space, an "expansion volume," 
where the cone of Cherenkov light expands. Fi- 
nally the photons enter a detector consisting of 
multi-wire proportional chambers filled with a 
mixture of TEA and CH 4 gases. No light focus- 
ing is used; this is called "proximity- focusing" £]. 
The scheme is shown in the upper left of Fig. 1 . 

There are 30 photon detectors around the cylin- 
der. They subtend the same azimuthal angle as 
the radiators, which are also segmented into 14 
sections along their length of the cylinder. The 
gap between the radiators and detectors, called 
the "expansion gap", is filled with pure N2 gas. 
The wire chamber design is shown in Fig. [l]. 

2.2. Radiators 

LiF was chosen over CaF2 or MgF2, both of 
which are transparent in the useful wavelength 
region, because of smaller chromatic error. Orig- 
inally all the radiators were planned to be 1 cm 
thick planar pieces. However, since the refractive 
index of LiF at 150 nm is 1.5, all the Cherenkov 
light from tracks normal to the LiF would be to- 
tally internally reflected as shown in Fig. || (top). 
We could have used these flat radiators, but we 
would have had to tilt them at about a 15° angle. 
Instead we developed radiators with striations in 
the top surface, called "sawtooth" radiators (?]], 
as shown in Fig. || (bottom). 

The overall radiator shape approximates a 
cylinder of radius 82 cm. Individual radiators 
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Figure 2. Sketch of a plane radiator (top) and a 
sawtooth radiator (bottom) . Light paths radiated 
from a charged track normal to each radiator are 
shown. 

are placed in 14 coaxial rings of 30 crystals each, 
centered around the beam line and symmetrically 
positioned about the interaction point. Only the 
inner four radiators are sawtooths.(The 30 crys- 
tals segments are parallel to the wire chambers.) 
Inter-crystal gaps are typically 50 /im. The crys- 
tals are attached to the exterior surface of a 1.5 
mm thick carbon fiber shell with a low outgassing 
epoxy. 

2.3. Photon Detectors 

Construction was carried out in a class 100 
clean room that was dehumidified below 35%. 
Granite tables were used that were flat over the 
entire surface of a photon detector module to bet- 
ter than 15 /im. 

The photon detectors have segmented cathode 
pads 7.5 mm (length) x 8.0 mm (width) etched 
onto G10 boards. The pad array was formed 
from four individual boards, with 24 x 80 pads, 
with the latter separated into two 40 pad sections 
with a 6 mm gap. Each board was individually 
flattened in an oven and then they were glued 
together longitudinally on a granite table where 
reinforcing G10 ribs were also glued on. There 
are 4 longitudinal ribs that have a box structure. 
Smaller cross ribs are placed every 12 cm for extra 
stiffening. The total length was 2.4 m. 

Wire planes were separately strung with 20 /im 
diameter gold plated tungsten with a 3% admix- 
ture of rhenium produced by LUMA; the wire 
pitch was 2.66 mm, for a total of 72 wires per 
chamber. The wires were placed on and subse- 



3 



quently glued to precision ceramic spacers 1 mm 
above the cathodes and 3.5 mm to the CaF2 win- 
dows, every 30 cm. We achieved a tolerance of 50 
/jm on the wire to cathode distance. The spacers 
had slots in the center for the glue bead. 

Eight 30 cm x 19 cm CaF2 windows were glued 
together in precision jigs lengthwise to form a 2.4 
m long window. Positive high voltage is applied 
to the anode wires, while - high voltage is put on 
100 fim wide silver traces deposited on the CaF 2 . 
The spacing between the traces is 2.5 mm. To 
maintain the ability of disconnecting any faulty 
part of a chamber, the wire HV is distributed 
independently to 3 groups of 24 wires and the 
windows are each powered separately. 

2.4. Electronics 

The position of Cherenkov photons is measured 
by sensing the induced charge on array of 7.5 mm 
x 8.0 mm cathode pads. Since the pulse height 
distribution from single photons is expected to 
be exponential at low to moderate gas gains ||, 
this requires the use of low noise electronics. Pad 
clusters in the detector can be formed from single 
Cherenkov photons, overlaps of more than one 
Cherenkov photon or charged tracks. In Fig. || 
we show the pulse height distribution for single 
photons, and charged tracks. We can distinguish 
somewhat between single photons hitting the pad 
array and two photons because of the pulse height 
shapes on adjacent pads. The charged tracks give 
very large pulse heights because they are travers- 
ing 4.5 mm of the CH4-TEA mixture. 

To have as low noise electronics as possible, a 
dedicated VLSI chip, called VA_RICH, based on 
a very successful chip developed for solid state 
applications, has been designed and produced 
for our application at IDE AS, Norway. We 
have fully characterized 3,600 64 channel chips, 
mounted on hybrid circuits. For moderate values 
of the input capacitance Cj n , the equivalent noise 
charge measured ENC is found to be about: 

ENC = 130e" + (9e"/pF) x C in . (1) 

Its dynamic range is between 450,000 and 900,000 
electrons, depending upon whether we choose a 
bias point for the output buffer suitable for signals 
of positive or negative polarity or we shift this 
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Figure 3. Pulse height distributions from pad 
clusters containing single photons (solid his- 
togram) and charged tracks (dashed histogram). 
The line shows a fit of photon data to an exponen- 
tial distribution. One ADC count corresponds to 
approximately 200 electrons. The charged track 
distribution is affected by electronic saturation. 
The wire voltage is +1.5 kV and the window volt- 
age -1.2 kV. 

bias point to have the maximum dynamic range 
for signals of a single polarity. 

In our readout scheme we group 10 chips in 
a single readout cell communicating with data 
boards located in VME crates just outside the 
detector cylinder. Chips in the same readout cell 
share the same cable, which routes control signals 
and bias voltages from the data boards and out- 
put signals to the data boards. Two VA_RICH 
chips are mounted using wire bonds on one hy- 
brid circuit that is attached via two miniature 
connectors to the back of the cathode board of 
the photon detector. 

The analog output of the VA_RICH is transmit- 
ted to the data boards as a differential current, 
transformed into a voltage by transimpedance 
amplifiers and digitized by a 12 bit differential 
ADC. These receivers are part of very complex 
data boards which perform several important 
analog and digital functions. Each board contains 
15 digitization circuits and three analog power 
supply sections providing the voltages and cur- 
rents to bias the chips, and calibration circuitry. 
The digital component of these boards contains a 
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sparsification circuit, an event buffer, memory to 
store the pedestal values, and the interface to the 
VME cpu. 

Coherent noise is present. We eliminate this 
by measuring the pulse heights on all the chan- 
nels and performing an average of the non-struck 
channels before the data sparsification step gj. 
The pedestal width (rms) changes from 3.6 to 
2.5 channels with and without this coherent noise 
subtraction, respectively. The total noise of the 
system then is ^500 electrons rms. 

2.5. Gas System 

The gas system supplies several distinct vol- 
umes. The systems must: supply CH4-TEA to 
30 separate chambers, supply super-clean N2 to 
the expansion gap, supply super-clean N2 to a 
sealed single volume surrounding all the cham- 
bers, called the electronics volume, since this is 
the region where the front-end hybrid boards are 
present. In addition we need to test CH4-TEA for 
the ability to detect photons and test the output 
N2 for purity. 

It is of primary importance that the gas system 
must NOT destroy any of the thin CaF2 win- 
dows. We use computerized pressure and flow 
sensors with PLC controllers. The gas system 
works great. N2 transparency is >99%. Nothing 
has been broken! 



nated by off-line pulse height thresholds adjusted 
to keep occupancy of each channel below 1%. Fi- 
nally we eliminate clusters of cathode pad hits 
that are extended along the anode wires, but are 
only 1-2 pads wide in the other direction. 

4.2. Cherenkov Images 

We show in Fig. || the hit pattern in the detec- 
tor for a Bhabha scattering event (e + e~ — > e + e~) 
for track entering the plane (left image) and saw- 
tooth (right image) radiators. The shapes of the 
Cherenkov "rings" are different in the two cases, 
resulting from refraction when leaving the LiF ra- 
diators. The hits in the centers of the images 
are produced by the electron passing the RICH 
MWPC. 



Figure 4. Hit patterns produced by the parti- 
cle passing through the plane (left) and sawtooth 
(right) radiators. 



3. OPERATING EXPERIENCE 

The detector has been in operation since 
September of 1999. All but -2% of the detector 
is functioning. We lost 1% due to the breaking of 
one wire after about one year of operation. We 
have also lost 2% of the electronics chips. 

4. OFF-LINE DATA ANALYSIS AND 
PHYSICS PERFORMANCE 

4.1. Noise filtering 

Coherent noise suppression and data sparsi- 
fication are performed on-line to eliminate the 
Gaussian part of the electric noise. A small non- 
Gaussian component of the coherent electric noise 
is eliminated off-line, by using an algorithm too 
complicated for use in the data board DSP. The 
incoherent part of non-Gaussian noise was elimi- 



4.3. Clustering of Hits 

The entire detector contains 230,400 cathode 
pads, which are segmented into 240 modules of 
24 x 48 pads separated by the mounting rails and 
anode wire spacers. We cluster pad hits in each 
module separately. Pad hits touching each other 
either by a side or a corner form a "connected 
region." Each charged track reconstructed in the 
CLEO-III tracking system |nj is projected onto 
the RICH MWPC and matched to the closest 
connected region. If the matching distance be- 
tween the track projection and the connected re- 
gion center is reasonably small and the total pulse 
height of the connected region sufficiently high 
we associate this group of hits with the track. 
Local pulse height maxima in the remaining con- 
nected regions, so called "bumps" are taken as 
seeds for Cherenkov photons. We allow the pulse 



height maxima to touch each other by corners if 
the pulse height in the two neighboring pads is 
small relative to both bump hits. Hits adjacent 
to the bumps on the sides are assigned to them 
in order of decreasing bump pulse height. 

To estimate the position of the photon conver- 
sion point we use the center-of-gravity method 
corrected for the bias towards the central pad. 
For many Cherenkov photons we are able to de- 
tect induced charge in only one pad. Since the 
pad dimensions are about 8x8 mm 2 , the posi- 
tion resolution in this case is 8/-\/l2 = 2.3 mm. 
For charged track intersections, which induce sig- 
nificant charge in many pads, the position reso- 
lution is 0.76 mm. The position resolution for 
Cherenkov photons which generate multiple pad 
hits is somewhere in between these two values. 
In any case, the photon position error is not a 
significant contribution to the Cherenkov angle 
resolution (see below). 

4.4. Corrections to the Track Direction 

The resolution of the CLEO-III tracking sys- 
tem is very good in the bending view (the mag- 
netic field is solenoidal in CLEO) [|l0| . The track 
position and inclination angle along the beam axis 
is measured less precisely, with the silicon vertex 
detector playing the dominant role. The rms of 
the observed RICH hit residual is 1.7 mm. Since 
the RICH hit position resolution is 0.76 mm as 
measured by the residual in the perpendicular di- 
rection, the RICH MWPC can clearly help in pin- 
ning down the track trajectory. This, in turn, 
improves Cherenkov resolution, especially for the 
flat radiators for which we observe only half of the 
Cherenkov image and thus are quite sensitive to 
the tracking error. The improvement is as much 
as 50% is some parts of the detector. 

4.5. Reconstruction of Cherenkov Angle 

Given the measured position of the Cherenkov 
photon conversion point in the RICH MWPC, the 
charged track direction and its intersection point 
with the LiF radiator, we calculate a Cherenkov 
angle for each photon-track combination (# 7 ). 
We use the formalism outlined by Ypsilantis and 
Seguinot [||, except that we adopt a numerical 
method to find the solution to the equation for 
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the photon direction, instead of simplifying it to 
a 4 th order polynomial, which would allow an ana- 
lytical solution, but at the expense of introducing 
an additional source of error. Furthermore, using 
our numerical method, we calculate derivatives of 
the Cherenkov angle with respect to the measured 
quantities which allows us to propagate the de- 
tector errors and the chromatic dispersion to ob- 
tain an expected Cherenkov photon resolution for 
each photon independently (<7g). This is useful, 
since the Cherenkov angle resolution varies sig- 
nificantly even within one Cherenkov image. We 
use these estimated errors when calculating parti- 
cle ID likelihoods and use them as to weight each 
photon when considering the average Cherenkov 
angle for a track. 

4.6. Performance on Bhabha Events 

We first view the physics performance on the 
simplest type of events, Bhabha events and then 
subsequently in hadronic events. The Cherenkov 
angle measured for each photon is shown in Fig. |^. 

We note that Bhabha events are very low 
multiplicity compared with our normal hadronic 
events. They have two charged tracks present 
while the hadronic events have an average 
charged multiplicity of approximately 10. In ad- 
dition, the hadronic events have on the average 
10 photons, mainly from ir° decays. All of these 
particles can interact in the calorimeter and the 
splash-back can hit the RICH photon detector. 

This Cherenkov angle spectrum for single pho- 
tons has an asymmetric tail and modest back- 
ground. It is fit with a line-shape similar to 
that used by for extracting photon signals from 
electromagnetic calorimeters |ll| . The functional 
form is 

P(0\6 e v P ,ag,a,n)= (2) 




Here is the measured angle, 9 exp is the "true" 
(or most likely) angle and ag is the angular res- 
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Figure 5. The measured minus expected 
Cherenkov angle for each photon detected in 
Bhabha events, (top) for plane radiators and 
(bottom) for sawtooth radiators. The curves are 
fits to special line shape function (see text), while 
the lines are fits to a background polynomial. 

olution. To use this formula, the parameter n is 
fixed to value of about 5. 

The data in Fig. || are fit using this signal shape 
plus a polynomial background function. We com- 
pare the results of these fits for the resolution pa- 
rameter ag as a function of radiator ring for data 
and Monte Carlo simulation in Fig. ^. The sin- 
gle photon resolution averaged over the detector 
solid angles are 14.7 mr for the flat radiator and 
12.2 mr for the sawtooth. 

The number of photons per track within a ±3ct 
of the expected Cherenkov angle for each photon 
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Figure 6. The values of the angular resolution 
for single photons for data compared with Monte 
Carlo simulation as a function of radiator ring. 
Sawtooth radiators are in rings 1 and 2, near the 
center of the detector. 

is shown in Fig. |?] and shown as a function of ra- 
diator ring in Fig. ||. Averaged over the detector, 
and subtracting the background we have a mean 
number of 10.6 photons with the flat radiators 
and 11.9 using the sawtooth radiators. 

The resolution per track is obtained by taking a 
slice within ±3<r of the expected Cherenkov angle 
for each photon and forming an average weighted 
by l/<Tg. These track angles are shown in Fig. ^. 

The rms spreads of these distributions are iden- 
tified as the track resolutions. We obtain 4.7 mr 
for the flat radiators and 3.6 mr for the sawtooth. 
The resolutions as a function of radiator ring are 
shown in Fig. [Io| 

The Cherenkov angular resolution is comprised 
of several different components. This include er- 
ror on the location of the photon emission point, 
the chromatic dispersion, the position error in the 
reconstruction of the detected photons and finally 
the error on determining the charged track's di- 
rection and position. These components are com- 
pared with the data in Fig. [□]. 

4.7. Performance on Hadronic Events 

To resolve overlaps between Cherenkov images 
for different tracks we find the most likely mass 
hypotheses. Photons that match the most hy- 
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Figure 7. The number of photons detected on 
Bhabha tracks (top) for plane radiators and (bot- 
tom) for sawtooth radiators. The dashed lines are 
predictions of the Monte Carlo simulation. 

pothcsis within ±3ct are then removed from con- 
sideration for the other tracks. To study the 
RICH performance in hadronic events we use 
D*+ ->■ 7r+L>°, D° -> K~% + events. The charge 
of the slow pion in the D* + decay is opposite to 
the kaon charge in subsequent D° decay. There- 
fore, the kaon and pion in the D° decay can be 
identified without use of the RICH detector. The 
effect of the small combinatorial background is 
eliminated by fitting the D° mass peak in the 
K + tt~ mass distribution to obtain the number of 
signal events for each momentum bin. 

Single-photon Cherenkov angle distributions 
obtained on such identified kaons with the mo- 
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Figure 8. The number of photons as a function 
of radiator ring. 



mentum above 0.7 GeV/c are plotted in Fig. 12 
Averaged over all radiators, the single-photon res- 
olution is 13.2 mr and 15.1 mr for sawtooth and 
flat radiators respectively. The background frac- 
tion within ±3cr of the expected value is 12.8% 
and 8.4%. The background-subtracted mean pho- 
ton yield is 11.8 and 9.6. Finally the per-track 
Cherenkov angle resolution is 3.7 mr and 4.9 mr. 

4.8. Particle ID Likelihoods 

For parts of the Cherenkov image for the saw- 
tooth radiator, and for tracks intersecting more 
than one radiator there are some optical path 
ambiguities that impact the Cherenkov angle cal- 
culations. In the previous section we bypassed 
this problem by selecting the optical path that 
produces the closest Cherenkov angle to the ex- 
pected one (6e Xp ) for the given particle hypothesis 
(h). There is some loss of information in this pro- 
cedure, therefore, we use the likelihood method 
to perform particle identification instead of the 
per-track average angle. The likelihood method 
weights each possible optical path by the opti- 
cal probability (P op t), which includes length of 
the radiation path and the refraction probabili- 
ties obtained by the inverse ray tracing method: 
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Figure 9. Track resolutions in Bhabha events, 
(left) for plane radiators and (right) for sawtooth 
radiators. 
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where, Lh is the likelihood for the particle hy- 
pothesis h (e, fi, 7T, K Or p), Pbackground IS the 

background probability approximated by a con- 
stant and P S ig n ai is the signal probability given 
by the line-shape defined previously. In principle, 
the likelihood could include all hits in the detec- 
tor. In practice, there is no point in inspecting 
hits which are far away from the regions where 
photons are expected for at least one of the con- 
sidered hypotheses (we use ±5cr cut-off). 
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Figure 10. Cherenkov angle resolutions per track 
as a function of radiator ring for Bhabha events. 

An arbitrary scale factor in the likelihood defi- 
nition cancels when we consider likelihood ratios 
for two different hypotheses. The likelihood con- 
veniently folds in information about values of the 
Cherenkov angles and the photon yield for each 
hypothesis. For well separated hypotheses (typi- 
cally at lower momenta) it is the photon yield that 
provides the discrimination. For hypotheses that 
produce Cherenkov images in the same area of 
the detector, the values of the Cherenkov angles 
do the job. Since our likelihood definition does 
not know about the radiation momentum thresh- 
old, the likelihood ratio method can be only used 
when both hypotheses are sufficiently above the 
thresholds. When one hypothesis is below the ra- 
diation threshold we use a value of the likelihood 
for the hypothesis above the threshold to perform 
the discrimination. 

The distribution of the 21n(L,r/ 'Lk), is ex- 
pected to behave as the difference \ 2 K — x\- This 
X 2 difference obtained for 1.0-1.5 GeV/c kaons 
and pions identified with the D* method is plot- 
ted in Fig. [H| Cuts at different values of this 
variable produce identification with different effi- 
ciency and fake rate. Pion fake rate for different 
values of kaon identification efficiency is plotted 
as a function of particle momentum in Fig. |lj. 

4.9. Conclusions 

We have successfully constructed and operated 
a large, complex RICH detector in a particle 
physics experiment for over three years. About 
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Figure 11. Different components of the 
Cherenkov angle resolutions per track as a func- 
tion of radiator ring for Bhabha events. The 
points are the data and the solid line is the sum 
of the predicted resolution from each of the com- 
ponents indicated on the figure. 

98% of the detector is operational (1% loss due 
to one broken wire and 1% due to electronic fail- 
ures) . 

The particle momenta for B meson decay prod- 
ucts seen by CLEO are less than 2.65 GeV/c. The 
detector provides excellent separation between pi- 
ons and kaons at and below this cutoff. Sep- 
aration between kaons and protons extends to 
even higher momentum, where it is used in charm 
baryon studies. Thus, the physics performance 
has met design criteria. 

CLEO currently is making an extensive study 
of Upsilon decays and proposes to study decays 
of charm mesons and charmonium decays (called 
CLEO-c |l2|]). For these measurements the beam 
energy will be lowered and the maximum particle 
momenta will be about 1.5 GeV/c. At these mo- 
mentum the particle identification fake rates are 
at the 1% level. 
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Figure 12. The measured minus expected 
Cherenkov angle for each photon detected in 
hadronic events, (top) for plane radiators and 
(bottom) for sawtooth radiators. The curves are 
fits to special line shape function (see text), while 
the lines are fits to a background polynomial. 
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Figure 13. Distribution of 2\n.(L^ / Lk) ~ x% ~ 
x\ for 1.0-1.5 GeV/c kaons (filled) and pions 
(open) identified with the D* method. 
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